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Abstract 

One possible feature of quantum gravity may be the violation of Lorentz invariance. 
In this paper we consider one particular manifestation of the violation of Lorentz 
invariance, namely modified dispersion relations for massive neutrinos. We show how 
such modified dispersion relations may affect atmospheric neutrino oscillations. We 
then consider how neutrino telescopes, such as ANTARES, may be able to place 
bounds on the magnitude of this type of Lorentz invariance violation. 
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1 Introduction 



The search continues for a theory which describes gravity on quantum scales, 
and provides us with falsifiable experimental predictions (for a review, see, for 
example, P). In recent years, the subject of quantum gravity phenomenology 
has rapidly grown (see [2] for a review), complementing theoretical work. 
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One particular feature of quantum gravity which has attracted much atten- 
tion in particle and astro-particle physics phenomenology is the breaking of 
CPT symmetry (see, for example, [3] for a review). CPT symmetry is antic- 
ipated to be broken in quantum gravity because the prerequisite conditions 
of quantum coherence, locality or Lorentz invariance may no longer hold. In 
previous work |^|f5|l6] we have studied the effect of quantum decoherence on at- 
mospheric and astrophysical neutrino oscillations, and in the present work we 
turn our attention to Lorentz invariance violation. We focus on atmospheric 
neutrino oscillations, having considered Lorentz invariance violating effects on 
astrophysical neutrinos in [6]. As in our previous work [4J, in this paper we 
take a phenomenological viewpoint. We examine, independent of any one par- 
ticular theory of quantum gravity, what we may expect to observe if Lorentz 
symmetry is broken or deformed. 

There is already a large body of work relating to Lorentz invariance violation 
and neutrinos, some of which is briefly reviewed in the following section. Fur- 
thermore, existing atmospheric neutrino data has been thoroughly analyzed 
for both Lorentz invariance violating effects and other new physics (see, for 
example, [TfSIU] ). Although the analysis of Super-Kamiokande and K2K data 
[7115] strongly disfavours Lorentz invariance violation as either the leading or 
sub-leading cause of atmospheric neutrino oscillations, an analysis of high- 
energy atmospheric neutrino oscillations, as would be observed at a neutrino 
telescope, would enable a wider range of Lorentz invariance violating models 
to be probed, as well as offering the possibility of further strengthening exist- 
ing bounds on Lorentz invariance violating parameters. This can be seen in 
the recent analysis of MACRO data [9] using both low (< 30 GeV) and high 
(> 130 GeV) energy neutrinos. 

In this paper, we consider the sensitivity of the ANTARES neutrino telescope 
[TU] to Lorentz invariance violating effects although we expect a similar anal- 
ysis to be applicable to other neutrino telescopes, such as ICECUBE [TT] . 
NESTOR [12] and NEMO [13]. Whilst neutrino telescopes are primarily de- 
signed to observe neutrinos from astrophysical and cosmological origins, at- 
mospheric neutrinos form the primary background to these sources and it is 
likely that they will provide the first physics results. 

The outline of this article is as follows. In section [2l we briefly discuss some 
of the theoretical consequences from Lorentz invariance violation, in partic- 
ular modified dispersion relations. We then model how modified dispersion 
relations alter atmospheric neutrino oscillations in section [31 We present the 
results of our simulations for the ANTARES neutrino telescope in section HJ 
where we first consider how the observed spectra of events may be altered by 
Lorentz invariance violating effects and then present the results of our sensi- 
tivity simulations. In section [5l we present our conclusions and compare the 
results obtained here with those from previous analyses of experimental data 
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[71l8ll9] . Unless otherwise stated, we use units in which c = h = 1. 



2 Quantum gravity violations of Lorentz invariance and the effect 
on neutrino oscillations 

In quantum gravity, space-time is expected to take on a foamy nature [l4j . 
Models of space-time foam have been built in many frameworks, including 
critical p^][T6] and non-critical string theory p2], and loop quantum gravity 
[18] . Interactions with this space-time foam may lead to the breaking of CPT 
symmetry, leading to the violation of Lorentz invariance [TS]. In addition, some 
theories of quantum gravity predict that there is a minimum length scale, of 
order the Planck length (10~^^ m) p3|l20] . If this is the case, then the existence 
of a fundamental length scale may also lead to Lorentz invariance violation 
(LV). 

To illustrate this, imagine if we were able to somehow measure this fundamen- 
tal length scale. We could measure it in our lab frame and probably find that it 
was of order the Planck length, 10~^^ m. If we now repeated the experiment in 
a different inertial frame (moving relative to the first frame), we would expect 
the length to be shorter, as predicted by the usual Lorentz contraction. How- 
ever, we are considering a fundamental length scale and, therefore, we would 
also expect to measure the same length in all inertial frames. There is, there- 
fore, either an error in our naive thought experiment or Lorentz symmetry, at 
least at the Planck scale, is broken. 

Therefore, there exist three possibilities for the fate of Lorentz invariance 
within theories of quantum gravity: 

(1) It may be that our naive thought experiment is wrong, or that a fun- 
damental length scale does not exist, and so Lorentz invariance remains 
intact at Planck scales; 

(2) It may be that, at the Planck scale, Lorentz invariance is broken and 
there is a preferred frame; 

(3) It may be that Lorentz invariance as we know it is deformed, so, for 
example, the Lorentz transformation involves a second quantity which is 
not observer dependent. So, perhaps the Planck length is invariant in an 
identical way to the speed of light. 

In the second and third options, one consequence of such breaking/ deforming 
of Lorentz invariance may be that quantities we think of as being Lorentz 
invariant, such as the standard dispersion relation, E"^ = p'^ + m^, may also 
have to be modified. 
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Lorentz invariance violation may manifest itself in many different ways and 
in many different experimental systems [2T|22] . Our focus here is on neutri- 
nos and modified dispersion relations. However, modifications of dispersion 
relations have also been considered for photons [23121], leading to dispersion 
relations of the form |24j : 



where uj is the frequency of the photon, k its momentum and M is the mass 
scale of the quantum gravitational theory. By considering photons which orig- 
inate in astrophysical objects such as gamma ray bursts and pulsars, lower 
bounds of ~ 10^^ GeV [21] were placed on M by considering the delay in 
arrival times. 

In the literature, the effects of LV on neutrino oscillations can be broken down 
into two distinct categories. In the first category are those models in which 
the only modification is to the neutrino oscillation length [6]|25y26f27j . which 
is the type of LV model we shall consider in this paper. 

In the second category are those models in which both the oscillation length 
and also the mixing angle are altered. In our simulations, we found that in- 
cluding LV effects on the mixing angle as well as the neutrino oscillation length 
made no difference to our results and so we do not study this type of model 
in detail. This type of LV, where both the mixing angle and oscillation length 
are altered, has been considered for the ICECUBE [28] and MINOS [29j ex- 
periments in addition to future neutrino factories [30] . 



3 Modified dispersion relations and neutrino oscillation 
probabilities 

Our interest in this paper is how LV effects, as described in the previous 
section, may alter the oscillations of atmospheric neutrinos. Therefore, in the 
analysis which follows, we restrict ourselves to a simple two neutrino system 
although it is possible to consider how these effects manifest themselves in a 
three neutrino system [6]. 

If we allow LV then this may lead to modifications of the dispersion relation 



where E is the energy of the neutrino, p the neutrino momentum, m is the 
mass eigenvalue of the neutrino and is the Planck energy. The function / 
contains all the novel LV effects. For our analysis here, we find it useful to 





E^=p^ + m^ + f{p,E,Ej,), 



(1) 
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parameterize this function so the dispersion relation ([T]) becomes [32] 



T7i2 2 I 2 I 2 

= p + m + rip 



E_ 

E„ 



(2) 



where r] and a are parameters of the Lorentz invariance violating theory. 



Using ([2]) and assuming the parameter rj has a dependence upon the mass 
eigenstate, then we find that the Hamiltonian in the mass basis, for a two 
neutrino system, may be written, to first order in t], as 



H 



2E ' 2 







2E ^ 2 



(3) 



where and rji are the mass and LV parameter of the ith eigenstate. Here we 
are considering only the case when the Hamiltonian is diagonal in the mass 
basis, as off-diagonal terms will modify the neutrino mixing angle as well as 
the oscillation length. 

The probability that a neutrino of flavour a oscillates to one of flavour b is 
then given by 

2 



-HiiJT* 
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where the f/'s are the components of the neutrino mixing matrix 

/ 



U 



V 



cos 9 sin 9 
— sin 9 cos 9 



with 9 being the mixing angle. The oscillation probability is therefore 
PWa Vb] = sin2 2^sin2 



AE 



AE'^ 
p 



where Am^ = m2 — mf and Arj = r]2 — rji. Replacing c and h and absorbing 
the Planck energy into the parameter Arj gives 



sin^ 29 sin^ 



1.27^^ + 1.27 X 109("+i)Ar7E"L 
E 



(4) 



where Am^ is measured in eV^, the neutrino energy E in GeV, the path length 
L in km, and we have set = a + 1. If we set Arj = in (jl]), then we recover 
the standard oscillation probability 



P\i^a ^b] = sin^ 29 sin^ 



1.27- 



.Am^L 
E 



(5) 
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whereas, if we set Am^ = 0, we find a probability in which neutrinos oscillate 
due to LV effects only: 




a 



Ub] = sin2 20sin2 \l.27 x 10^^''+^^ ArjE^'L . 



In ([3]), we have assumed that the LV parameter rj has a dependence upon the 
mass eigenstate. If this is not the case and the parameter rj is universal to 
all mass eigenstates, then, even if the theory violates Lorentz invariance, the 
neutrino oscillation probability is unchanged. 

In this paper, we shall consider three models with differing energy depen- 
dences, namely: 

(1) Model LVl - with n = 1: here the LV parameter Arj in (jlj) is dimension- 
less; 

(2) Model LV2 - with n = 2: then Ar] has units of eV"^; 

(3) Model LV3 - with ri = 3: in which case Arj has units eV~^. 

It may be, in a phenomenology of quantum gravity, that n is not an integer 
but we only consider integer values here (the results for non-integer n can 
be found by interpolating between the two bounding integer values). Models 
with n = —1, and 1 have been examined using Super-Kamiokande [7] data, 
and more complicated models with n = or 1 have been analyzed using 
both Super-Kamiokande and K2K data [8]. We do not consider n < 1 here 
as n = models have been severely constrained by these previous analyses, 
and, since in that case the LV modifications of the oscillation probability are 
independent of the neutrino energy, we do not expect to be able to improve 
on their bounds by using higher-energy atmospheric neutrinos. Therefore we 
only consider values of n which are strictly positive, and include n = 1 for 
comparison with previous work [7115] . From our work on quantum decoherence 
effects on atmospheric neutrino oscillations one might anticipate that we 
will be able to place more stringent bounds on LV effects as n increases. 

In the next section, we shall be concerned with finding upper bounds on the 
LV parameters. However, we are able to extract crude order of magnitude 
estimates for these parameters from the oscillation probability (jl]). The LV 
effects only become significant when 



The peak in the atmospheric flux is of the order of 1 GeV, which is the 
typical energy scale for experiments such as Super-Kamiokande [33], whilst 
neutrino telescopes, such as ANTARES, are sensitive to much higher energy 
neutrinos. Table [1] shows a comparison of the typical bounds on Arj for low 
and high energy neutrino experiments. Here, we have taken the peak in the 
atmospheric flux, E ^ 1 GeV as the value of the neutrino energy for low 
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F — 1 PpV 


F _ 1 fin 


1 


2.3 X 10"2i 


2.3 X 10-25 


2 


2.3 X 10-3° eV-i 


2.3 X 10-36 eV-i 


3 


2.3 X 10-3^ eV-2 


2.3 X 10"^^ eV-2 



Table 1 

Table showing the expected upper bounds on Ary, the LV parameter for various 
values of n, for atmospheric neutrinos with low energies, E ^ 1 GeV, and higher 
energies, E ~ 100 GeV. 

energy experiments and ~ 100 GeV as a typical value for neutrinos to 
be detected by a neutrino telescope such as ANTARES. As a comparison, 
the "most conservative" bound for n = 1 from Super-Kamiokande data f7] is 
6 X 10"^'^ at the 90% confidence level. The authors of [7J comment that the 
strongest bounds are coming from upgoing muon events, which peak at about 
100 GeV, so their results are in line with naive expectations. The bounds of 
[7] are improved by a factor of about 8 in the combined analysis of Super- 
Kamiokande and K2K data [8] . We would hope that the much greater number 
of such events that will be observed at a neutrino telescope would enable these 
bounds to be improved. 

To see how Lorentz invariance violation affects atmospheric neutrino oscilla- 
tions, it is useful to compare the oscillation probabilities. Whilst the neutrino 
oscillation probabilities themselves are not observables, they are useful for 
gaining insight into the underlying physics. Firstly, for comparison, we con- 
sider the standard oscillation probability which is given by (jj]). In all the 
probabilities presented here, we take [M] 

Am^ = 2.6 X 10-3 eV^; sin^ 26 = 1. 

Fig. [1] shows the oscillation probability plotted as a function of the path length, 
L, for neutrino energy 200 GeV and as a function of the neutrino energy, E, 
for a fixed path length, L = 10'^ km (which is comparable with the diameter 
of the Earth, and used in the oscillation probabihty plots in [9j). 

Figs. [21, [3] and H] show the oscillation probability including LV effects (jl]) for 
n = 1, 2 and 3 respectively. To construct these plots, we used the upper 
bounds found in section H] and shown in table [2l In each case, we found 
that for low energy neutrinos, the effects of LV on the neutrino oscillation 
probability are negligible. However, for higher energy neutrinos, this is not 
the case. If we compare the plots of probability against path length for high 
energy neutrinos, we see that for standard oscillations, almost no oscillations 
occur, however, if we include the LV effects, the oscillation probability be- 
comes significant at much shorter path lengths. Indeed, as n increases, the 
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Fig. 1. Standard atmospheric neutrino oscillation probability ([5]). The left plot shows 
the oscillation probability as a function of path length L (measured in km) , for fixed 
neutrino energy E equal to 200 GeV. The right plot show the oscillation probability 
as a function of energy for path length L = 10^ km. 
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Fig. 2. Atmospheric neutrino oscillation probabilities for n = 1 (a = 0). The left 
plot shows the oscillation probability as a function of path length for neutrinos with 
an energy of 200 GeV whilst the right plot shows the oscillation probability as a 
function of energy for path length L = 10^ km. 
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Fig. 3. As Fig. [2] but with n = 2 (a = 1). 

probability becomes significant at smaller path lengths. When we consider the 
plots in which the probability is a function of energy, we see that for stan- 
dard oscillations, the probability tends towards zero, where as, in the plots 
which include LV effects, the probability continues to oscillate with the period 
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Fig. 4. As Fig. [2] but with n = 3 {a = 2). 

of oscillations decreasing for increasing values of n. Corresponding plots of 
atmospheric neutrino survival probabilities can be found in P] for the mass 
and velocity eigenstate mixing LV model. Their plots have some features in 
common with ours, in particular the existence of oscillations in the probabil- 
ity at high energies, when the survival probability in the standard oscillation 
scenario has converged to unity. 

In |1] , we studied in depth the effect of quantum decoherence on atmospheric 
neutrino oscillations. In that case, there were also marked differences in the 
oscillation probability, particularly for high energies and long path lengths. 
However, the probabilities with LV are rather different from those when we 
have quantum decoherence. In the decoherence case, the oscillations in the 
probability are damped (the rate of damping depending on a number of fac- 
tors), but here we find that the oscillations in the probability continue for ever. 
While these two possible sources of new physics give very different properties 
of the oscillation probability, the key test is whether these differences can be 
observed experimentally. We turn to this issue in the next section. 



4 ANTARES sensitivity to LV in atmospheric neutrino oscillations 

The ANTARES neutrino telescope [lU] is sensitive to high energy neutrinos 
above the energy threshold of 10 GeV. These neutrinos do not form the pri- 
mary source of neutrinos for ANTARES, instead they form the background 
to astrophysical and cosmological neutrinos. 

4-1 Simulations 

Our simulation strategy follows that presented in [4j, so we only briefly outline 
the key features. Further details are documented in [1], and in the analysis of 
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standard atmospheric neutrino oscillations [35], on which our simulations are 
based (with suitable modifications for the inclusion of LV). 

Atmospheric neutrino events were generated using a Monte-Carlo (MC) pro- 
duction assuming a spectrum of E""^ for neutrinos in the energy range 10 GeV < 
E < 100 TeV. The angular distribution was assumed to be isotropic. A to- 
tal of twenty five years of data were simulated in order that errors from the 
Monte-Carlo statistics could be ignored and event weights are used to adapt 
the Monte-Carlo fiux to a real atmospheric neutrino fiux. The Bartol theoreti- 
cal flux [36] was used here although we do not anticipate the results presented 
will change using a different flux, however, there may be a signiflcant change 
to our results if a different spectral index is assumed. 

The simulations presented here are based on three years' worth of data taking 
with ANTARES, corresponding to around 10000 atmospheric neutrino events. 
All errors are purely statistical with simple Gaussian errors assumed. We are 
able to produce spectra in either E, E/L or L and, by altering the oscillation 
probability in the neutrino oscillation code, we are able to include LV effects. 
The path length L is related to the zenith angle (not to be confused with 
the mixing angle 6) as follows, provided cos-i? ^ 0: 

L ^ 2i?cos^9, 

where R = 6378 km is the radius of the Earth, and i) = corresponds to a 
neutrino travelling vertically upwards. We are therefore interested in typical 
path lengths of the order L 10^ km. ANTARES is particularly sensitive to 
the zenith angle and hence the path length L. In the sensitivity analysis, we 
therefore used the spectra in E/L. However, the spectra in L alone also have 
some interesting features, which we discuss below. 

In order to produce sensitivity regions, we use a technique which compares 
the foi' oscillations with LV, or LV effects alone, with that for the no- 
oscillation hypothesis. We leave the total normalization as a free parameter 
so that our results are not affected by the normalization of the atmospheric 
neutrino flux. Sensitivity regions at both the 90 and 99 percent confldence 
level are produced. 

Although our analysis is independent of the normalization of the atmospheric 
neutrino flux, in order to produce spectra based on the number of observed 
events, some normalization is necessary. For standard oscillations with no LV 
effects, bins at high E/L can be used to normalize the flux as standard oscil- 
lations become negligible here. However, since all our LV models are propor- 
tional to positive powers of the neutrino energy, LV effects are signiflcant at 
high energies but negligible at low energies, and thus low E/L. In this case, 
we use the low E/L bins to normalize the flux. 
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In our simulations, we consider cases where Am^ is zero, so neutrino oscil- 
lations arise due to LV effects only, and where Am^ is non-zero, so that the 
LV effects modify standard neutrino oscillations. In the case when Am^ is 
non-zero, we also allow this parameter to vary in order to check whether the 
experimental point of best fit lies inside our sensitivity regions. In this case, 
we have three parameters, Am^, Ar/ and sin^26'. We therefore obtain a three 
dimensional sensitivity region at 90 and 99% confidence levels. In order to 
consider the details of the sensitivity region, we plot the boundaries of the 
projection of the surface bounding this volume onto the coordinate planes, 
thus obtaining three distinct plots. 



4-2 Spectra 

In addition to data based on the total number of observed events, the spectra 
of these events, typically in E/L, is also important. For example, spectral 
analysis has been performed for Super-Kamiokande [7f51l5B] . KamLAND [H7] 
and K2K ^8i5B] . In addition to being used in the construction of sensitivity 
regions, the spectra themselves are also important for bounding the size of 
possible LV effects and possibly ruling effects out due to the observed shape of 
the spectra. As explained above, we used the spectra in E/L in the sensitivity 
analysis, so we focus on those spectra (apart from Fig. 

4.2.1 Model LVl 

We begin by examining the spectra obtained when we consider the possibility 
that neutrino oscillations occur as a consequence of LV only {Am? = in Eq. 
(jl])), when the LV effects are proportional to the neutrino energy, n = 1. Fig. 
[5] shows the expected spectra of events as a function of E/ cos-i? and cos-i?. 
The solid line corresponds to the MC simulation with standard oscillations 
whereas the dotted/dashed lines are the MC simulations for LV only induced 
oscillations. In each case, we have also plotted simulated data points, based 
on the standard oscillation picture with no LV effects. On the right-hand-side, 
we show the ratios of the number of events compared to no oscillations. 

The spectra show that, for the larger values of the LV parameter, there is a 
definite suppression in the number of observed events by a factor of approx- 
imately two. If, however, we consider the case in which the LV parameter 
takes the smaller value, it is more difficult to distinguish oscillations that arise 
as a consequence of LV effects from those arising in the standard oscillation 
picture. If we consider the ratios, we note that there are some features which 
would allow us to differentiate between the two phenomenologies. In the top 
right frame, we note that at low energies, the LV model does not result in an 
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Fig. 5. Spectra of events for atmospheric neutrino oscillations due to LV effects 
only, when the LV effects are proportional to the neutrino energy (left panel) and 
their ratios to the events without oscillations (right panel) as functions oi E/ cos 
(upper plots) and cosi? (bottom plots). The solid line represents the MC simulation 
of standard oscillations with the dotted/dashed lines showing the spectra for os- 
cillations from LV only. The simulated data points correspond to three-years' data 
taking, assuming standard oscillations only with no LV. 

oscillation minimum whereas the standard oscillation model does. Secondly, 
although the difference is small, the bottom right frame shows a deviation 
from the standard picture at high values of the zenith angle corresponding 
to large path lengths. We will not consider further spectra just in the path 
length cos??, but note that for all the models we study and all the energy 
dependences of the LV parameters, there are slight deviations in the shape 
of the spectra in cosd compared with the standard oscillation case. In the 
analysis of Super-Kamiokande data (7|l8| it is the spectra in the path-length L 
which are used to bound the magnitudes of LV parameters, as the inclusion 
of LV effects again alters the shape of the spectra. 

For the remainder of this section, we focus on the ratios of the numbers of 
events in each case compared with the expected number of events if there were 
no oscillations. Fig. [6] shows the ratios in the case of standard oscillations plus 
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Fig. 6. Ratio of the number of events for standard oscillations modified by LV effects 
proportional to the neutrino energy, compared with no oscillations. The solid line 
represents the MC simulation of standard oscillations with the dotted/dashed lines 
showing the spectra for standard oscillations modified by LV. The simulated data 
points assume only standard oscillations, with no LV effects. 

LV as a function oi E / cosd. Again, we note that large values of the LV pa- 
rameters result in a large suppression of the number of expected events. Also, 
for very large values of the LV parameter, we see a considerable flattening in 
the ratio of the number of events compared with no oscillations. The situation 
is much trickier, however, for very small values of the LV parameter. In this 
case, the number of events expected is not suppressed and we also see an oscil- 
lation minimum. Therefore, this makes it very difficult to distinguish between 
standard oscillations and those which are modified to include LV effects, if the 
LV parameter is of the order of 10~^^ or smaller. 




4.2.2 Models LV2 and LV3 

We first consider the spectra when the LV parameters are proportional to the 
square of the neutrino energy. Fig. [7] shows the ratio of number of events with 
the LV effects only to those expected if there were no oscillations. For large 
values of the LV parameters, we note that the spectra are very different from 
standard neutrino oscillations. Firstly, as in the case before, the number of 
events is significantly suppressed for very large values of the LV parameter. For 
standard neutrino oscillations, we observe an oscillation minimum, however, 
for large values of the LV parameter, the ratio increases and so we see a 
peak. Even for smaller values of the LV parameter, the ratio is fiat and so by 
considering the spectra it should be possible to distinguish between these two 
models. 
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Fig. 7. Ratio of the number of events for oscillations due to LV effects only, when 
the LV effects are proportional to the neutrino energy squared, compared with no 
oscillations. The solid line represents the MC simulation of standard oscillations 
with the dotted/dashed lines showing the spectra for oscillations from LV only. 
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Fig. 8. As Fig. [7] but the dotted/dashed lines show the spectra for standard oscilla- 
tions plus LV. 

Fig. [8] shows the spectra in the case where the LV effects modify standard 
oscillations, then, as in subsection 14.2.11 before, it becomes much more diffi- 
cult to untangle any LV effects from those of standard oscillations for small 
LV parameter. For large values of the LV parameter, the only distinguishing 
feature is that the number of events we expect to see is somewhat reduced. In 
addition, for large values of the LV parameter, the curve is much flatter than 
we would expect from standard oscillations. 

We find qualitatively similar results when the LV effects are proportional to 
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Fig. 9. Sensitivity contours at 90 and 99% confidence levels for LV effects only (no 
standard oscillations) for model LVl, when the LV effects are proportional to the 
neutrino energy. 

the cube of the neutrino energy, both in the case where oscillations are due 
to LV effects only, and when LV effects modify standard neutrino oscillations. 
For both models LV2 and LV3, the shape of the spectra in the zenith angle 
cos only are slightly different to that shown in for model LVl in Fig. [51 The 
ratios of the number of events compared with no oscillations are flatter than 
in the bottom right-hand plot in Fig. 0, particularly for model LV3. We do not 
use these spectra for the following sensitivity analysis, but it may be possible, 
with real data, to use those spectra as well as those for numbers of events in 
terms oi E / cosd. 



4-3 Sensitivity regions 

In the previous section, we saw how LV effects can modify the number of events 
seen in the detector and how they would modify the spectra. We now turn to 
the discussion of the ANTARES sensitivity regions found from our numerical 
simulations. As in the case of quantum decoherence [4J, we are particularly 
interested in finding upper bounds for the LV parameters. 



4.3.1 Models LVl and LV2 

In our first model, the LV effects are proportional to the neutrino energy. We 
firstly consider the case in which Am^ = 0, so that there are no standard 
oscillations, only those which arise as a consequence of the LV effects. Fig. M 
shows the sensitivity curves in this case. It is clear that the smallest values of 
the LV parameter may be probed when sin^ 29 is close to one, which is the 
region of parameter space in which we are most interested. We are able to 
place an upper bound on Ar] of 8.2 x 10~^^. 
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Fig. 10. Sensitivity contours for model LVl at 90 and 99% confidence level for stan- 
dard oscillations plus LV effects proportional to the neutrino energy. The triangle 
denotes the best-fit value of Am^ with no LV effects [HI] . 

We also considered the case of non-zero Am^, when standard atmospheric 
neutrino oscillations are modified by LV effects. Fig. [10] shows the sensitivity 
contours in this case. Since we have varied three parameters (sin^(26'), Am^ 
and Ar/), we obtain a sensitivity volume. Projecting the boundary surface 
of this volume onto the co-ordinate planes gives us the sensitivity contours 
plotted in Fig. [TOl As we see from Fig. [TOl the experimental point of best 
fit of Am^ [34J, denoted by the triangle, lies within the sensitivity region. 
However, the region extends down to very small values of Am^ and so we find 
that our results include the previous case where neutrino oscillations arise as 
a consequence of LV effects only. The top right frame in Fig. [TU] shows the 
sensitivity contour for the Arj as a function of the mixing angle, 9. Again, 
At] = is in this region, and so standard oscillations with no LV effects are 
also included, with an upper limit shown in Table [2] This is the same order of 
magnitude as the upper bound from Super-Kamiokande and K2K data [7]l8] 
so it seems unlikely that ANTARES will be able to significantly improve on 
current data. 



We have performed a similar analysis for the model LV2, when the LV effects 
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n 


Lxr] [e\ ) 

( Am'^ — (Tl 
\L^lli — U ) 


L\ri [e\ ) 


Planck scaling 
{e\ ) 


1 


8.2 X 10~25 


2.9 X 10^24 


X 


2 


1.0 X 10-35 


2.9 X 10^35 


10-28 


3 


X 


6.9 X 10""^^ 





Table 2 

Table showing the upper bounds on Aij contained with the sensitivity regions for 
various values of n. We also give the naive expected values of Ar] from suppression 
by appropriate powers of the Planck mass. The X indicates that we were unable to 
place a bound on this parameter. 

are proportional to the neutrino energy squared. The sensitivity regions for 
the cases where oscillations are due to LV effects only, and where standard 
oscillations are modified by LV effects, are very similar in shape to those shown 
in Figs. [9] and [10] respectively. In particular, the experimental best-fit point 
shown by a triangle in Fig. [10] again lies within the sensitivity region when 
standard oscillations are modified by LV effects. In this case, the sensitivity 
region again extends to both Am^ = and Arj = 0. The upper bound on Arj 
from this sensitivity region is given in Table [2] 

4.3.2 Model LV3 

The final model we consider has n = 3. If we set Am^ = 0, so that oscillations 
are due to LV effects only, then we find that we are unable to derive any 
meaningful sensitivity contours, so we do not consider this possibility further. 
For Am^ ^ 0, the sensitivity contours are shown in Fig. [11] The sensitivity 
contours in this case have a rather simpler shape than in models LVl and 
LV2. Again the experimental point of best fit is included in the sensitivity 
region but in this case, Am^ = is not included, so it seems that we are not 
sensitive to the case when oscillations arise from LV effects only for this energy 
dependence. This is in direct agreement with the first part of our analysis as we 
were unable to find any meaningful results for this model when we considered 
LV effects only. We note however that Arj = is contained in the sensitivity 
region. The upper bound on At] in this sensitivity region is given in Table [2] 



5 Discussion and conclusions 

We have shown that neutrino telescopes, such as ANTARES, will be able to 
place stringent bounds on LV parameters which modify the dispersion rela- 
tion for massive neutrinos. Table [2] shows the upper bounds for this model 
for the various energy dependences we have studied. We found that if the LV 
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Fig. 11. Sensitivity contours for model LV3 at 90 and 99 percent confidence level 
for standard oscillations plus LV effects proportional to the neutrino energy cubed. 
The triangle denotes the best-fit value of Am'^ with no LV effects |34] . 

effects were proportional to either the neutrino energy or the neutrino en- 
ergy squared, then we are sensitive to oscillations induced from these effects 
only, with no contributions from standard oscillations. However, for the case 
of n = 3, then we cannot measure the effects of LV only in atmospheric neu- 
trino oscillations. For all values of n studied, we have been able to produce 
sensitivity curves for standard atmospheric neutrino oscillations modified by 
subdominant LV effects. We should stress that although our simulations are 
based on the ANTARES neutrino telescope, we expect that other neutrino 
telescopes, such as AMANDA and ICECUBE, will also be able to probe these 
models. 

From a naturalness perspective, we would expect any LV effects to be sup- 
pressed by the Planck mass, 10^^ GeV. Since our LV parameters Arj contain 
powers of the Planck mass already, this gives Ar] ~ l0-28{n-i) oyi-n n>2. 
These values are also given in table [2]for comparison. We note that for the case 
n = 2, the bounds obtained from our simulations surpass the expected value 
of the LV parameter by seven orders of magnitude and so ANTARES may be 
able to observe these effects or rule out this model. For the case n = 3, how- 



18 



ever, the bound from our simulations is nine orders of magnitude larger than 
expected and so we cannot rule out this model using atmospheric neutrinos. In 
this case better sensitivity can be obtained using higher-energy, astrophysical 
neutrinos [H]. The situation is somewhat more complicated if n = 1 as this 
corresponds to a = in ([2]) and so we have no indication as to the value of Arj 
except to say that it should be small. In all three cases, we can see from table 
|2] that our estimates based on the values of the LV parameters at which LV 
effects become comparable to standard neutrino oscillations are very similar 
to the upper bounds we find from our sensitivity curves. 



Our results for the case n = 1 are comparable with the model investigated 
in jTlIB] with Super-Kamiokande and K2K data. We find that we are not able 
to significantly improve on their results; however an analysis of atmospheric 
neutrino data from ANTARES should be able to confirm their upper bounds. 
Although in a rather different LV model, our results are also comparable in 
magnitude to those found by the MACRO collaboration P|. In that paper, 
the argument of the sine term in the oscillation probability is identical to 
that considered here and they found an upper bound on Ar] in the range 
10~^^ < Af] < 10~^^, similar to that found here and shown in table [2l 
Our results are also comparable to those found in [27] using data from the 
AMANDA experiment. There, an upper bound of 2.1 x 10^^'' was placed on 
the LV parameter. 



In [1], we studied in depth the spectra of numbers of events for standard at- 
mospheric neutrino oscillations with quantum decoherence. Those spectra had 
some of the same key features as the ones presented here in Figs. [MH in par- 
ticular, a suppression in the expected number of atmospheric neutrino events 
observed by ANTARES. The shape of the spectra do have some differences, 
but it remains to be seen whether they will really be distinguishable exper- 
imentally. At this point, in the absence of real data, one can only speculate 
on what might be observed, but one can foresee two possibilities. Firstly, it 
may be that ANTARES (or any other neutrino telescope) does not observe 
any distinguishable difference in numbers of events compared with the stan- 
dard atmospheric neutrino oscillation model. We argued in [4J that such a null 
observation would be able to place very stringent bounds on quantum deco- 
herence effects, and we draw the same conclusion here from our event spectra. 
Secondly, it may be the case that ANTARES or an equivalent neutrino tele- 
scope does indeed find evidence for a suppression in the number of events 
compared with the standard picture. In this situation it is likely to be very 
challenging experimentally to disentangle the possible source of new physics. 
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